Generation of High Quality Laser Accelerated Ion Beams 
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In order to achieve a high quality, i. e. monoergetic, intense ion beam, we propose the use of a 
double layer target. The first layer, at the target front, consists of high-Z atoms, while the second 
(rear) layer is a thin coating of low-Z atoms. The high quality proton beams from the double 
layer target, irradiated by an ultra-intense laser pulse, are demonstrated with three dimensional 
Particle-in-Cell simulations. 
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The high efficiency of ion acceleration recently ob- 
served in the interaction of petawatt laser pulses Q with 
solid targets || has lead to important applications such 
as the use of laser produced fast ion beam injection to 
conventional accelerators (see Ref. |J), the fast ignition 
of thermonuclear targets, as discussed in Refs. Q], and 
hadrontherapy in oncology |5j . Laser accelerated protons 
have been used for proton imaging of small scale objects 
in laser produced plasmas with high time resolution ||. 
The process of ion acceleration has been studied in de- 
tail with multi-dimensional Particle-in-Cell (PIC) simu- 
lations 0. In the recent experimental results presented 
in Refs. |!|, electron energy in a range of hundreds of 
MeV was observed, while the proton energy was about 
tens of MeV with a number of fast protons from 10 12 
to 10 13 per pulse and with the transformation of 12% 
of the laser energy into fast ion energy. In Refs. it 
was shown with PIC simulations that, by optimizing the 
laser-target parameters, it becomes possible to accelerate 
protons up to several hundreds MeV with a number of 
fast ions approximately 10 13 particles per pulse. 
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FIG. 1: Electric field near the target is shown as the 3D vector 
field, the length of each vector corresponds to the magnitude 
of electric field. Vectors with \eE x / (m e u)c)\ < 1 are violet. 
The laser pulse is presented by isosurfaces of the transverse 
component E z corresponding to dimensionless values ±10. 
The time unit is the laser period 2-k/uj. 
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FIG. 2: Plasma species. The envelopes of heavy ions (thick 
red shell) and light ions (thin blue shell) are shown. Electron 
density is shown as a 'green gas' using ray tracing technique. 



The typical energy spectrum of laser accelerated par- 
ticles observed both in the experiments and in the 
computer simulations can be approximated by a quasi- 
thermal distribution with a cut-off at a maximum en- 
ergy £ max . The effective temperature T, that may be at- 
tributed to the fast ion beams, is much less than the max- 
imal energy. On the other hand, almost all above men- 
tioned applications require high quality proton beams, 
i.e., beams with sufficiently small energy spread A£/£. 
For example, for the hadron therapy it is highly desirable 
to have a proton beam with A£/£ < 2% in order to pro- 
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FIG. 3: The proton (blue) and heavy ion (red) energy spectra 
at t = 80 • 2-k/u. 

vide the conditions for a high irradiation dose being deliv- 
ered to the tumor saving neighboring healthy tissues |^| . 
In the concept of Fast Ignition with laser accelerated ions 
presented in Refs. 0, the proton beam was assumed to 
be quasi-monoenergetic. An analysis carried out in Ref. 

has shown that the Fast Ignition with a quasi-thermal 
beam of fast protons requires several times larger energy 
than that with a monoenergetic beam. Similarly, in the 
case of the injector (see Ref. H), a high-quality beam is 
needed in order to inject the charged particles into the 
optimal accelerating phase efficiently. Thus we see that 
the generation of high quality beams is a key problem for 
many applications. However energy spectra of the laser 
accelerated ions at present are rather far from required 
ones. 

In this Letter we show with three dimensional (3D) 
PIC simulations that such a required beam of laser ac- 
celerated ions can be obtained using a double layer target 
(see also Refs. ||). Multi layer targets have been used for 
a long time in order to increase the efficiency of the laser 
energy conversion into plasma and fast particle kinetic 
energy (see Refs. in In contrast to the configu- 

rations discussed previously, we propose to use a double 
layer target to produce fast proton beams with controlled 
quality. 

In the proposed scheme the target is made of two lay- 
ers. The first layer consists of high-Z atoms (atomic mass 
Ttii ) , while the second layer is a very thin coating of low- 
Z atoms (atomic mass ro a ). Such the target can be a 
metal foil coated with a thin hydrogeneous film. An 



ultra-intense laser pulse is incident on the first layer, and 
so we say that the first layer is at the target front, while 
the second layer is at the rear side of the target. We use 
the term 'longitudinal' for the direction of propagation 
of the laser pulse, and the term 'transverse' for perpen- 
dicular directions. 

When an ultra- intense laser pulse irradiates the target, 
heavy atoms are partly ionized and electrons are expelled 
from the foil. A quasi-static electric field is generated 
due to charge separation. The first layer of heavy ions 
(the foil) should be sufficiently thick so as to produce a 
large enough quasi-static electric field, and, at the same 
time, it should be sufficiently thin to be able to produce 
strong electric field at the rear side. Such an electric 
field has opposite signs on the two sides of the target, 
and vanishes inside the target and at some finite distance 
from it. The number of low-Z ions in the second layer (the 
coating) should be sufficiently small not to produce any 
significant effect on the electric field. The quasi-static 
electric field accelerates both high-Z ions (with average 
charge eZ{) and low-Z ions (with average charge eZ a ). 
If the ratio rriiZ a / (m a Zi) is sufficiently large, light ions 
are accelerated much more efficiently than heavy ones. 
The coated thin layer can be then accelerated forward. 
It detaches from the foil, and moves as a whole in the 
longitudinal direction. The light ions within a small solid 
angle have a quasi-monoenergetic energy spectrum. The 
thinner the coating the narrower energy spectrum of light 
ions is. 

The important requirement is that the transverse size 
of the coating must be smaller than the laser waist since 
an inhomogeneity in the laser pulse causes nonuniform 
accelerating electric field and thus a degradation of beam 
quality, as seen in experiments presented in Refs. 
where the exposed targets had a thin proton layer on 
their surface. The laser pulse inhomogeneity results in 
an additional energy spread of the ion beam as seen in 
the experiments. The effect of the finite waist of the laser 
pulse leads also to an undesirable defocusing of the fast 
ion beam. In order to compensate for this effect and 
to focus the ion beam, we can use properly deformed 
targets, as suggested in Refs. [0, [lT| . 

In a target, whose transverse size is much greater than 
the laser focal spot, quasi-static electric field can be af- 
fected by background 'cold' electrons from the periphery. 
To increase the quality and life-time of the quasi-static 
electric field, we should use targets with transverse size 
less than the laser waist. 

In order to estimate the typical energy gain of fast 
ions, we assume that the main portion of the free elec- 
trons produced by ionization in the irradiated region of 
the foil are expelled. In this case the generated elec- 
tric field near the positively charged layer is equal to 
Eq = 2TmoeZil. Here no is the ion density and eZi is 
an average ion electric charge in the foil, and I is the 
foil thickness. The region of strong electric field has a 
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transverse size of the order of the diameter 2R± of the 
focal spot. Thus the longitudinal size of the region where 
the electric field remains essentially one-dimensional is of 
the order of R± and the typical energy of an ion with 
charge eZ a accelerated by this electric field can be es- 
timated as £ max = 2nn Z a Zie 2 lR±. We have assumed 
that the electron energy in the laser field is well above 
the ion energy and larger than the energy required for the 
electrons to leave the irradiated region. The maximum 
electron energy in the electromagnetic wave is given by 
£ e = m e a 2 /2, where a = eE / (m e cuj) is the dimensionless 
amplitude of the laser pulse. From this condition we can 
find the required value of the laser pulse intensity and its 
power. 

Consider a double layer target in the form of a prolate 
ellipsoid coated at the rear side with a very thin proton 
layer. We assume that ellipsoid semiaxes are R± and 1/2. 
We can estimate the energy spectrum of protons, using 
formulae for the electric field of the electrically charged 
prolate ellipsoid Jlj|. Let the cc-axis be in the longitu- 
dinal direction, originating from the target center. On 
this axis the x-component of the electric field is given by 
E x (x) = (E /3)R 2 J(R 2 ± - (l/2) 2 +x 2 ). The distribution 
function of the fast protons f(x,v,t) obeys the kinetic 
equation d t f + vd x f + (eE(x)/m p ) d v f — 0. If particles 
trajectories do not intersect or self-intersect (for our ellip- 
soidal target this it true), we have f(x,v,t) — fo(xo,Vo), 
where fo(xo,vo) is the distribution function at the ini- 
tial time t = 0. The number of particles per unit vol- 
ume in phase space, dxdv, is equal to dn = fdxdv = 
fvdvdt = f d£ dt/rrip. We assume that at t = all 
particles are at rest and that their spatial distribution 
is given by no(xo) which corresponds to the distribu- 
tion function fo(xo,Vo) = no(xo)5(vo), with S(vq) the 
Dirac delta function. For the applications discussed in 
this paper we are interested in the particle distribu- 
tion function integrated over time. Time integration 
of the distribution fvdvdt gives the energy spectrum of 
the beam N(£)d£ = {no{xo)/m p )\dt/dv\ v=Vo d£. Here 
the Lagrange coordinate of the particle xq and the Ja- 
cobian \dt/dv\ v=v are functions of the particle energy 
£. The Lagrange coordinate dependence on the energy 

= x o(£) is given implicitly by the integral of the par- 
ticle motion: £(x,xq) = £o + e[ip(x) — ip(xo)], where ip(x) 
is the electrostatic potential. In the case under consid- 
eration, we have £q = and x = oo. The Jacobian 
\dt/dv\ v=v is equal to the inverse of the particle accel- 
eration at t = 0, i.e. \dt/dv\ v=VQ = l/\eE x (xo)\. On the 
other hand the function \dt/dv\ v=Vo is equal to \dxo/d£ \. 
Hence, we obtain the expression for the energy spectrum 
in the form 

N(£)d£ = n (x )d£ /\d£ /dx \ Xo=Xo{£ y (1) 

We notice that the expression for the energy spectrum 
follows from the general condition of particle flux conti- 
nuity in the phase space. 



As we can see, in the vicinity of the target the elec- 
tric field on the axis is homogeneous and equals to 
E x (l/2) — 2TrnoZiel/3. Therefore, the form of the en- 
ergy spectrum (|l|) is determined from the distribution of 
the proton density n (x = tp~ 1 (£/e)). We see that in 
general a highly monoenergetic proton beam can be ob- 
tained when the function no(^o) is a strongly localized 
function, i.e. when the thickness of the proton layer Axq 
is sufficiently small, and then A£/£ as Axo/R±. 

The longitudinal emittance of the beam is defined as 
the product of the energy spread A£ and the time length 
At of the beam: en = A£At. Using the expressions 
obtained above we find the longitudinal omittance of the 
accelerated proton beam e| = (Axo/R±) 2 y/m p £R 2 L /2. 
For £ = lOOMeV, R± = 5^m and Ax /Rj_ = 6 ■ 1CT 3 we 
obtain e| w 1.3 ■ 10~ 4 MeV ps which is about 150 times 
smaller than the emittance observed in the experiments 
with non-optimized targets in Ref. p2| . 

Near the axis, the radial component of the electric field 
depends linearly on the radius r = \/y 2 + z 2 : E r (r) = 
(EqI i '6i?jjr . We find that the particle trajectory is de- 
scribed by r = ro exp(y/kx), where Tq is the initial radial 
coordinate of the particle and k = l/Rj_. From this ex- 
pression, for l/R± <C 1 we find the transverse emittance 
e± = nroAd of the fast proton beam: e± = Tvroy/l/R±. 
Here ro is the transverse size of the proton layer and 
A9 = cxp(\/fci?jj — 1 is the beam divergence angle. For 
ro ~ 2.5/im, = 5/im and I « 0.5yum the transverse 
omittance is of the order of 2.5 mm mrad. 

For the hadron therapy the particle flux must be ap- 
proximately 10 10 to 5 • 10 10 protons per second ||. If 
the laser pulse is focused onto a spot with diameter 
10/im, CH layer with average density lg/cm 3 and thick- 
ness 0.03/^m provides 10 11 fast protons per pulse. 

In order to take into account the numerous nonlinear 
and kinetic effects as well as to extend our consideration 
to a multidimensional geometry, we performed numerical 
simulations of the proton acceleration during the inter- 
action of a short, high power laser pulse with a two-layer 
target. We used the three dimensional massively parallel 
and fully vectorized code REMP (Relativistic Electro- 
Magnetic Particle-mesh code) |l4| . In these simulations 
the largest number of grid cells was 2560 x 1024 x 1024 
and the number of quasiparticles was up to 820- 10 6 . The 
boundary conditions for the particles and for the fields 
are periodic in the transverse directions and correspond 
to absorption at the ends of the computation box along 
the x axis. The simulations were performed on 64 pro- 
cessors of the vector supercomputer NEC SX-5 at CMC, 
Osaka University. 

Here we present the results of these simulations. The 
size of the simulation box is 80A x 32A x 32A. A lin- 
early polarized laser pulse with dimensionless amplitude 
a = 30 propagates along the x-axis. The pulse size is 
15A x 12A x 12A. The pulse has a trapezoidal shape 
(growth-plateau-decrease), with 3A — 2A — 10A in the x- 
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direction, and 1A — 10A — 1A in the y- and z-directions. 
The plasma consists of three species: electrons, protons 
with m p /m e — 1836, and heavy ions (gold with Z,- L = +2) 
with mi/m e Zi = 195.4 x 1836/2. 

The first layer (gold) is placed at x = 5.5A. It has the 
form of a disk with diameter 10A and thickness 0.5A. The 
second layer (protons) also has the form of a disk with 
diameter 5A and thickness 0.03A, and is placed at the 
rear of the first layer, at x — 6A. The electron density in 
the heavy ion layer corresponds to the ratio uj pe /oj = 3.0 
between the plasma and the laser frequencies, for the 
proton layer it corresponds to w pe /w = 0.53. The number 
of electrons in the first layer is approximately 2000 times 
larger than that in the proton layer. 

The simulation results are shown in Figs. |l| to ||, where 
the coordinates are measured in wave lengths of the laser 
light and the time in laser periods. In Fig. [j] we present 
the electric field inside the computation box, to show the 
shape of the transmitted laser pulse and the longitudinal 
electric field accelerating the protons. The electric field is 
shown as a three-dimensional vector field; it is localized 
in the vicinity of the first layer (the layer of heavy ions) of 
the target and can be described as an electrostatic field 
from a positively charged disk. The transmitted laser 
pulse is presented by the isosurfaces of constant value of 
the z-component of the electric field. In Fig. || we show 
the densities of plasma species inside the computation 
box. We see that the proton layer moves along the x-axis 
and that the distance between the two layers increases. 
The heavy ion layer expands due to Coulomb explosion 
and tends to become rounded. Part of the electrons is 
blown off by the laser pulse, while the rest forms a hot 
cloud around the target. We notice that for the simula- 
tion parameters the electrons do not expelled from the 
region irradiated by the laser light completely. Even if 
only a portion of the electrons is accelerated and heated 
by the laser pulse, the induced quasi-static electric field 
appears to be strong enough to accelerate the protons up 
65 MeV, as seen in Fig. g presenting the energy spectra 
of the protons and heavy ions. The energy per nucleon 
acquired by the heavy ions is approximately 380 times 
smaller than the proton energy. The heavy ions have a 
wide energy spectrum while the protons form a quasi- 
mono-energetic bunch with A£/£ < 5%. We emphasize 
that the thinner the proton layer (or low-Z coating) the 
better the beam quality is. The proton beam remains lo- 
calized in space for a long time due to the bunching effect 
of the linearly decreasing electric field in the acceleration 
direction. 

In conclusion, the use of the multi-layer targets with 
different shapes and compositions opens up new oppor- 
tunities for controlling and optimizing the parameters of 
the fast proton (ion) beam, such as the energy spectrum, 
the number of particles per bunch, the beam focusing 
and the size of the region where the beam deposits its 



energy. 
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